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Abstract 

We compute the survival probabihty of large rapidity gaps (LRG) in a QCD based eikonal 

model with a dynamical gluon mass, where this dynamical infrared mass scale represents the onset 
of nonperturbative contributions to the diffractive hadron-hadron scattering. Since rapidity gaps 
can occur in the case of Higgs boson production via fusion of electroweak bosons, we focus on 
WW — > H fusion processes and show that the resulting (jS'p) decreases with the increase of the 
energy of the incoming hadrons, in line with the available experimental data for LRG. We obtain 
(IS'P) = 27.6 ± 7.8 % (18.2 ± 7.0 %) at Tevatron (CERN-LHC) energy for a dynamical gluon mass 
mg = 400 MeV. 
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I. INTRODUCTION 



The study of the survival probabihty (|5'p) of large rapidity gaps (LRG) is currently 
a subject of intense theoretical and experimental interest. Its importance lies in the fact 
that systematic analyses of LRG open the possibility of extracting New Physics from hard 
diffractive processes. On the theoretical side, its significance is due to the reliance of the 
(1 5*1^) calculation on subtle QCD methods. 

Rapidity gaps are defined as regions of angular phase space devoid of particles jjj, 
Ol , and at high-energy hadron-hadron collisions it has been suggested that their observation 



may serve as a signature for production of 



3iggs bosons and other colour-singlet systems via 



fusion of electroweak bosons |l[ 0, 0, flSQB. However, as pointed out by Bjorken QQ, 
we are not able to distinguish between the theoretically calculated rate of a large rapidity 
gap, Fgap, and the actual measured rate, fgap. These rates are related by the proportionality 
relation fgap = {\S\'^) Fgap, where is the so called "survival probability" of a large 

rapidity gap: it gives the probability of survival of LRG and not the probability for the 
production and survival of LRG, which is the quantity actually measured. More specifically, 
the survival factor gives the probability of a large rapidity gap not to be filled by 

debris originated from the rescattering of spectator partons {{\Sspec\'^)), or from the emission 
of bremsstrahlung gluons from partons {{\Shre,n\'^)). Hence we can write 

im = {\S,reU^y=\y,-y2m{\S,pUs)\'), (1) 

where yi and y2 are the rapidities of the jets (produced with large transverse momenta), 
= \yi ~ y2\ ^ l, and s is the square of the total center-of-mass (CM) energy. The 
factor {\Sbrem{^y)\'^) , which depends on the value of the LRG as well as on the kinematics 
of each specific process, can be calculated using perturbative QCD jssl. On the other hand, 
the survival factor (|S'specP) cannot be calculated perturbatively: it takes account of soft 
rescatterings in both the initial and final state interactions. To calculate (IS'^pecP) we need 
to obtain the probability of the two incoming hadrons do not interact inelastically, i.e. the 
probability of large rapidity gaps do not be populated by additional hadrons produced from 
the rescattering of spectator degrees of freedom. As discussed by Gotsman, Levin and Maor 
some time ago, this task is a difficult one since partons at long distances contribute to such 
a calculation However, the survival probability can be properly defined in the impact 
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parameter space 

,^_ J£bF{b,s)P{b,s) 

^ " Jd%F{b,s) ' ^ ^ 

where b is the impact parameter, F{b, s) is a factor related to the overlap of the parton 
densities in the colliding hadrons in the transverse impact plane, and P{b,s) is the probability 
that the two incoming hadrons have not undergone a inelastic scattering at the parton 
level. The above definition sets up the stage for carrying out a systematic computation 
of (IS'specP), since the probability P(6,s) can be easily obtained, for example, from the 



QCD-inspired eikonal approach Ul, 112; lid; llJl • ^^^^ picture the probability that neither 
hadron is broken up in a collision at impact parameter b is given by P{b, s) = e-^x/C'.^)^ where 
Xi{b,s) is the imaginary part of the eikonal function. In QCD-inspired eikonal models the 
increase of the total cross sections is associated with semihard scatterings of partons in 
the hadrons, and the high energy dependence of the cross sections is driven mainly by 
gluon-gluon scattering processes. Nevertheless, the gluon-gluon subprocess cross section is 
potentially divergent at small transferred momenta, and the usual procedure to regulate this 
behaviour is the introduction of a purely ad hoc parameter separati ng t he perturbative from 
the non-perturbative QCD region, like an infrared mass scale jlll . Il^ , or a cut-off at low 
transverse momentum px 

Recently we introduced a QCD-based eikonal model where the ad hoc infrared mass scale 
was substituted by a dynamical gluon mass one Q]. One of the central advantages of the 
model is that it gives a precise physical meaning for the quoted infrared scale. Furthermore, 
since the behaviour of the running coupling constant is constrained by the value of dynamical 
gluon mass jl^^^; the model also has a smaller number of parameters than similar QCD 
models. 

In this letter we perform a detailed computation of the survival probability in pp and pp 
channels in the framework of the QCD eikonal model with a gluon dynamical mass. We 
are concerned with the calculation of {ISspecl"^), the probability that LRG survive the soft 
rescattering of spectator partons, which we shall denote henceforth simply as (l^p). In the 
next section we introduce the QCD eikonal model and address the question of calculating 
the survival factor (jS'p) from processes of Higgs boson production through W fusion. The 
results are presented in the Sec. Ill, where we provide a systematic study of (|5'p) and its 
sensitivity to the infrared mass scale. The conclusions are drawn in Sec. IV. 



II. LARGE RAPIDITY GAPS AND THE DYNAMICAL GLUON MASS 



It has been suggested that in hadron-hadron colhsions the production of Higgs bosons 
can occur by means of the fusion of gluons or electroweak bosons . The two main 

Higgs production mechanisms are the gluon fusion gg H and the W fusion WW H. 
In the gluon fusion process each hadron emits a gluon (a colour octet), and the Higgs boson 
is coupled to each one through a fermion loop. Since the Higgs boson couples to fermions 
according to their masses, its production cross section via gluon fusion is dominated by top 
quark loops. In gluon fusion the hadrons remnants must exchange colour with each other 
in order to become singlet states again. 

On the other hand, the W fusion process has a different colour flow structure: since the 
incoming hadrons (as well as the W bosons) are colour singlets, when they emit a W boson 
they remain as singlet states. Thus no colour is exchanged between the scattered partons, 
and the two outgoing hadron remnant states are expected to be separated by a central 
rapidity gap. As indicated in the previous section, the survival probability of these rapidity 
gaps can be naturally defined in the impact parameter representation. In this formalism the 
inclusive differential Higgs boson production cross section via W fusion is given by 



d(yprod 



aww^HW{b;ij,w), (3) 



cm 

where W{b; fiw) is the overlap function at impact parameter space of the W bosons. This 
function represents the effective density of the overlapping W boson distributions in the 
colliding hadrons. The cross section for producing the Higgs boson and having a large 
rapidity gap is given by 

= aww^H W{b; ^w) P{h, s), (4) 

where P(&, s) is the probability that the two initial hadrons have not undergone a inelastic 
scattering at the parton level. In QCD-inspired eikonal models this probability is given 
by P(6, s) = e-^^^^^^''^), where the imaginary part Xii^i^) of the eikonal function receives 
contributions of parton-parton interactions. Therefore, the factor P(6, s) suppresses the 
contribution to the Higgs boson cross section where the two initial hadrons overlap and 
there is soft rescatterings of the spectator partons. From the expressions © and Q we 



4 



can write down the survival factor ([S'p) for Higgs production via W fusion: 



d%Wib;fiw)e-^'''^'^'\ (5) 



where we have used the normahzation condition J d'^bW{b; fiw) = 1- Iii this letter we shall 
compute the probability factor P{b, s) = e"^-^^*^^'*-* using a recently developed QCD eikonal 
model, where the onset of the dominance of gluons in the interaction of high-energy hadrons 
is managed by the dynamical gluon mass scale 1^. The model, henceforth referred to 
as DGM model, satisfies analyticity and unitarity constraints. The latter is automatically 
satisfied in the eikonal representation, where the total cross section, the ratio p of the real to 
the imaginary part of the forward scattering amplitude, and the differential elastic scattering 
cross section are given by 

atotis) =471 f bdb[l- e~^r(b''^ cosx^ib, s)], (6) 
Im{^/W6[l-e*x{M)]}' 



and 



^(.,t) = i- jbdb\l-e^^^''^^]Mqb{ 



(8) 



respectively, where s is the square of the total CM energy, Jo{x) is the Bessel function of the 
first kind, and s) = Xr{^i s) + ixA^i s) is the (complex) eikonal function. In the DGM 
model the eikonal function is written as a combination of an even and odd eikonal terms 
related by crossing symmetry. In terms of the proton-proton {pp) and antiproton-proton 
{pp) scatterings, this combination reads x%ip, s) = x^ib, s) ± x i^, s). The even eikonal is 
written as the sum of gluon-gluon, quark-gluon, and quark-quark contributions: 

X^{b,s) = Xqq{b,s) + Xq9{b,s) + X9g{b,s) 

= i[agq{s)W{b; /x^g) + aqg{s)W{b; pgg) + agg{s)W{b; /x^g)], (9) 

where W{b;fi) is the overlap function at impact parameter space and o"ij(s) are the ele- 
mentary subprocess cross sections of colliding quarks and gluons = q,g). The overlap 



5 



function, normalized so that J d'^bW{b; n) = 1, is associated with the Fourier transform of 
a dipole form factor, 

W{b;fi) = ^ifib)'K,if^b), (10) 

where K^lx) is the modified Bessel function of second kind. The odd eikonal x~{b,s), that 
accounts for the difference between pp and pp channels, is parametrized as 

X-(6,s) = C-S^e-/^iy(6;/x-), (11) 

where is the dynamical gluon mass and the parameters and fi^ are constants to be 
fitted. The factor S is defined as 



(12) 



with the dynamical coupling constant set at its frozen infrared value. The origin of the 
dynamical gluon mass and the frozen coupling constant can be traced back to the early work 
of Cornwall [ISj, and the formal expressions of these quantities can be seen in Ref. Iq . 

The eikonal functions Xqq{b, s) and Xggi^, s), needed to describe the lower-energy forward 
data, are simply parametrized with terms dictated by the Regge phenomenology: 

X,,(6,s)=zSA^l^(6;/iJ, (13) 



Xqgib,s) = 



A' + B'\n 



s 

rrig 



W{b;,/JI^), (14) 



where A, A' , B', figg and fitting parameters. The gluon-gluon eikonal contribution, 

that dominates at high energy and determines the asymptotic behaviour of the total cross 
sections, is written as Xggip, s) = cr^^^ {s)W {b; figg), where 

^r^^) = C'f drF,,{r)a^^\s). (15) 

J im^ I s 

Here Fgg^r) is the convoluted structure function for pair gg, cTgg^is) is the subprocess 
cross section and C is a fitting parameter. In the above expression it is introduced the 
energy threshold s > 4m^ for the final state gluons, assuming that these are screened gluons 
jisl ]. The structure function Fgg{T) is given by 

Fggir) = [g0 g]iT) = 9{^) 9 (^) , (16) 



where g{x) is the gluon distribution function, adopted as 



9[x) 



(17) 



The correct analyticity properties 
_5g-*7r/2 throughout Eqs. (US)), ^ 



yDPT 
99 



fx 



is calculated 



where J = 1 + e and Ng = — e)(5 — e)...(l — e) 

of the model amplitudes is ensured by substituting s - 
and (US)). 

In the expression fllSj) the gluon-gluon subprocess cross section 
using a procedure dictated by the dynamical perturbation theory (DPT) [19|: amplitudes 
that do not vanish to all orders of perturbation theory are given by their free-field values, 
whereas amplitudes that vanish in all orders in perturbation theory as oc exp {—l/g"^) {g is 
the coupling constant) are retained at lowest order. In our case this means that the effects 
of the dynamical gluon mass in the propagators and vertices are retained, and the sum of 
polarizations is performed for massless (free-field) gluons. As a result, since the dynamical 
masses go to zero at large momenta, the elementary cross sections of perturbative QCD in 
the high-energy limit are recovered. Other details of the calculation can be seen in Ref . . 

According to the expression (jSl), the final step in order to calculate (IS'p) is to determine 
the overlap function W{b] fiw) of the electroweak bosons. It is worth mentioning that the 
survival factor (ISp) depends on the nature of the colour- singlet exchange which generates 
the gap as well as on the distributions of partons inside the proton in impact parameter 



space 



20 



21 



22l |23[ . We simply assume that the distribution of W bosons in impact 



parameter space in the hadron is the same as for the quarks. In this way, we can finally 
write down a phenomenologically useful expression to the survival factor (l^p): 



{\S\^) = 2n / bdbW{b;fXgg)e 



-2X1 {b,s) 



(18) 



In the above expression the inverse size (in impact parameter) /i^g is the same as in 
the expression ()13|) . Its value, as well as the value of the remaining fitting parameters, is 
determined from global fits to pp and pp forward scattering data, as discussed in the next 
section. 



III. RESULTS 

In order to determine the exponential damping factor e~'^^^^^''^^ of Eq. ()18p and produce 
a consistent estimate of (l^'p), we carry out global fits to the elastic differential scattering 
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cross section for pp at ^/s = 1.8 TeV and to all high-energy forward pp and pp scattering 
data above y/s = 15 GeV. This energy threshold is the same one used in the estimate of 
(1 5*1^) through the analysis of m and pp scattering carried out by Block and Halzen using a 
previous QCD-inspired model The forward data sets include the total cross section (atot) 
and the ratio of the real to imaginary part of the forward scattering amplitude (p) . We use 
the data sets compiled and analyzed by the Particle Data Group [2.4], with the statistic and 
systematic errors added in quadrature. The input values of the rrig have been chosen to lie 
in the interval [350, 650] MeV, as suggested by the value m 



g — iidH,-ioo 



previous analysis of the pp and pp channels via the DGM model 



00+?^° MeV obtained in a 



IJj. This input dynamical 



gluon mass range is also supported by recent studies on the •yp photoproduction and the 
hadronic 77 total cross sections j^], and on the behaviour of the gluon distribution function 
at small x 26|. In all the fits performed in this letter we use a fitting procedure, assuming 
an interval ~ Xmm corresponding, in the case of normal errors, to the projection of the 

hypersurface containing 90% of probability. In the case of the DGM model (8 fitting 
parameters) this corresponds to the interval ~ Xmin = 13.36. 

The x^/DOF values obtained in the global fits are relatively low, as show in Table I. 
These results (for 168 degrees of freedom) indicate the excellence of the fits and show that 
the DGM model naturally accommodates all the data sets used in the fitting procedure. 
In Table I we have included the values of the ^w{= fJ'qq) parameter, which determines the 
spatial distribution of the W bosons at the impact parameter b. We can observe a small 
dependence of their values on the dynamical gluon mass: the greater the input scale irig, 
the smaller the inverse size fiw 

The sensitivity of the survival probability (for pp collisions) to the gluon dynamical 

mass is shown in Figure 1 for some CM energies. We note a slow increase of their values with 
the dynamical gluon mass and a fast decrease with the CM energies. As shown in Figure 
2, where we have also plotted the exponential damping factor, this behaviour is related to 
the energy dependence of the imaginary part of the eikonal: Xi{b, s) grows with the energy 
and hence suppresses the integral ()18|1 . In Table II we list our results for the survival factor 
(l^p) for some values of the proton-proton energy, and compare with other calculations in 
the literature, where {\S\^)DGAn and {\S\^) DGM2 denote the results obtained by setting the 
mass infrared scale at nig = 400 and 600 MeV, respectively. The last value is the same 
one adopted by Block and Halzen with respect to the ad hoc mass scale mo We see 
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that the DGM results are systematically larger than the Block-Halzen ones (denoted by 
bh), but the large statistical errors resulting from our choice for the confidence region 
of the parameters (CL=90 %) indicate a reasonable compatibility between the {\S\'^)dgmi 
and 

BH results at higher energies. 
The {\S\'^)£)GM2 results at y/s = 1.8 and 16 TeV are in line with the {\S\'^)glmi ones, 
obtained by Gotsman and collaborators using a Regge pole model j^^. In their approach 
;hey use an eikonalized version of the Donnachie-Landshoff model in order to satisfy unitarity 
3|. The authors argue that their relatively large values for (jSp) can be reduced by an 
appropriate change in some parameters included in their Gaussian approximation for F{b, s) 
and P{b, s) factors. We hasten to emphasise that the {\S\'^)r)GMi results have been obtained 
using the preferred statistical value of the dynamical gluon mass for pp and pp scattering, 
namely rrig = 400 GeV ^ . In this case we believe that an eventual change in the parameters 
of the GLMl model may reduce their results in such a way to be compatible with the DGMl 
ones. 

The {\S\'^)kmr results have been obtained by Khoze and collaborators using a two-channel 
eikonal model which embodies pion-loop insertions in the Pomeron trajectory, diffractive dis- 
sociation, and rescattering effects [2^. The authors have calculated the survival probability 
(ISp) in single, central and double diffractive processes at several energies, assuming that 
the spatial distribution in the parameter space is controlled by the slope b of the Pomeron- 
proton vertex. We show the {\S\'^)kmr results for double diffractive processes with 2b = 5.5 
GeV^, which corresponds to the slope of the electromagnetic proton form factor. These re- 
sults are compatible with the DGM ones, in particular with the results taking into account 
ra, ^ 400 MeV. the optimum value for the dyuamica. g.uon mass m aud » diffractive 
scattering |15| . 



IV. CONCLUSIONS 

In this letter we have calculated the survival probability (jS'p) of large rapidity gaps by 
means of an eikonal QCD model with a dynamical gluon mass. Since rapidity gaps can 
occur from production of Higgs boson via fusion of electroweak bosons, we have focused 
on WW H fusion processes. The eikonal function have been determined by fitting pp 
and pp accelerator scattering data. Owing to the quality of the global fits, the DGM model 
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allows us to describe successfully the pp differential cross section at y/s = 1.8 TeV, as well as 
the forward scattering quantities a^J^f^ and pPf^Pf^ in excellent agreement with the available 
experimental data. These results show that the DGM model is well suited for the prediction 
of the survival probability of LRG at higher energies, in particular for (ISp) one at the 
CERN-LHC energy. 

In Table II we list our results for and notice that their values decrease with the 

increase of the energy of the incomin g pr otons. This behaviour, in line with results for LRG 



dijet production at the Tevatron |30|, is a direct consequence of the energy dependence 



of the imaginary part of the eikonal, that grows with the energy. A strong dependence of 
(1 5*1^) on the dynamical gluon mass rrig emerges from our calculations, as shown in Figure 
1. This scale dependence arises as follows: the dynamical gluon mass affects strongly the 
behaviour of the gluon-gluon subprocess cross section o"^^^(s), which dominates at high 
energy and determines the asymptotic behaviour of the pp and pp total cross section. Hence 
the procedure consisting of global fits to diffractive pp and pp data inorder to determine rUg 



is well justified, and the value rrig = 400^^oo MeV obtained in Ref. [1^ via the DGM model 
is a suitable one. 

Our estimates for the survival probability of large rapidity gaps using a QCD based 
eikonal model with a dynamical gluon mass are, within the errors, compatible with estimates 
obtained using other eikonal models. In particular, our estimates are close to the ones 
obtained by Khoze et al. using a two-channel model, and to the ones obtained by Block 
and Halzen using a similar QCD-inspired approach. Owing to the interval [300,750] MeV 
inferred from the optimal value of mg discussed above, there is room for smaller values of the 
survival factor in DGM model. For example, a mass scale mg ~ 300 MeV gives a survival 
factor (IS'p) ~ 15.3 % at LHC, very close to the central value obtained via the KMR model. 

However, we call attention to the fact that all these estimates are model dependent, 
despite their apparent agreement. For example, the {\S\'^)kmr results for other values of 
26 and for central and single diffractive processes do not agree with ours [2^. The same is 
expected in the {\S\^)bh results for other choices of the mass scale uiq. 

In summary, there is a strong dependence in the size of the survival probabilities and 
in their energy dependence on specific models for the rise of total hadronic cross section. 
More specifically, the survival factor {\S\'^) depends on the dynamics of the whole diffractive 
part of the scattering matrix as well as the nature of the colour-singlet exchange which 
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generates the gap. From the experimental viewpoint, it is known that the survival factor 
in the case of Higgs production via WW H fusion processes can be monitored by 
observing the closely related central production of a Z boson with the same jet and rapidity 
gap configuration |3^. More recently, this idea has been developed further by considering 



the decays of both Higgs and Z bosons into bb pairs [33] • This allows to gauge Higgs weak 
boson fusion production at the LHC and permits to observe a light Higgs boson via its 
dominant H —>-bb decay mode in addiction to the usually discussed rr and WW* channels. 
This option would permit to reduce the theoretical uncertainty in the rate of Higgs central 
production events with rapidity gaps. 

The success of the QCD-based eikonal model with a dynamical gluon mass in reproducing 
diffractive scattering data, over a large energy range, shows that such a model provides a 
reliable estimate of the survival probability (l^p) as a function of energy in the case of pp 
and pp channels. The study of the survival factor (l^p) is interesting in its own right since 
they enables us to increase our understanding of some features of hadronic interactions, and 
may provide an useful tool to probe physics beyond the Standard Model. 
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TABLE I: The fiw parameter as a function of the dynamical gluon mass nig. The x^/DOF values 
resulting from the global fits are obtained for 168 degrees of freedom. 



mg [GeV] 




liw [GeV] 




XVDOF 


350 




0.8308±0.1394 




1.043 


400 




0.8091±0.1410 




1.022 


450 




0.7848±0.1411 




1.010 


500 




0.7823±0.1392 




1.009 


550 




0.7227±0.1356 




1.000 


600 




0.7254±0.1333 




1.001 


650 




0.7025±0.1305 




0.999 


TABLE II: The survival probability {\S\'^) (in %) for pp collisions in different models. 


^fs [GeV] 


{\S\^)dgmi 


{\S\'')dgm2 {\S?)bh 


{\S\'^)glmi 


{\S?)kmr 


63 


45.4±8.4 


50.9±9.3 37.5±0.9 






546 


34.2±8.1 


39.4±8.9 26.8±0.5 




26.0 


630 


33.4±8.1 
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FIG. 1: The survival probability {\S\'^) (central values) as a function of the dynamical gluon mass 
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FIG. 2: The imaginary part xi(b,s) of the eikonal and the exponential factor e~'^^'^^'^'^ for pp 
colUsions as a function of the impact parameter b, where = 1.8 TeV and ^/sJJ = 14 TeV. The 
dynamical gluon mass scale was set to mg = 400 MeV. 
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